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A finite element method to model dynamic structural 
systems undergoing large rotations is presented. The dynamic 
systems are composed of rigid joint bodies and flexible beam 
elements. The configurations of these systems are subject to 
change due to the relative motion in the joints among 
interconnected elastic beams. A body fixed reference is 
defined for each joint body to describe the joint body's 
displacements. Using the finite element method and the 
kinematic relations between each flexible element and its 
corotational reference, the total displacement field of an 
element, which contains gross rigid as well as elastic 
effects, can be derived in terms of the translational and 
rotational displacements of the two end nodes. If one end of 
an element is hinged to a joint body, the joint body's 
displacements and the hinge degree of freedom at the end are 
used to represent the nodal displacements. This results in a 
highly coupled system of differential equations written in 
terms of hinge degrees of freedom as well as the rotational 
and translational displacements of joint bodies and element 
nodes . 
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Introduction 


1 . 


In many space applications, such as deployment of large 
space lattice structures, maneuvers of space crane system, 
slewing of space antenna, and so on, the structural systems 
are composed of articulated elastic links which may undergo 
large angular motion relative to one another. The analytical 
simulation of the operation of such systems is critical to 
their successful design since the size and gravity suspension 
effects limit the quality of ground test while the flight 
experiments are still very expensive. This analytical 
simulation technique is referred to as multibody dynamic 
analysis . 

During the past decade, many efforts (Shabana and 
Wehage, 1984), (Singh and eial, 1985), (Yoo and Haug, 1986) 
have been devoted to the 3-D flexible multibody dynamic 
analysis. Because of the different characteristics of the 
mechanisms and/or . structural systems, various assumptions and 
theories were adopted, and different approaches were employed 
to formulate the equations of motion governing the systems' 
behavior. Since the configuration of an elastic component in 
a system is defined by infinite spatial degrees of freedom; 
theoretically, a set of partial differential equations of 
time and spatial variables are required to represent the 
system's dynamic equations. In order to take advantage of 
modern computational facilities, component mode synthesis 
techniques, (Shabana and Wehage, 1983), (Singh and ei al, 1985), 

and (Yoo and Haug, 1985) , have been used to approximate the 



dynamic model resulting in a set of numerically acceptable 
ordinary differential equations with finite spatial degrees 
of freedom. In most of the previous studies, the gross rigid 
motion of an elastic body is often defined by the 
translational and rotational displacements of a body 
reference frame which may be either fixed to or floating 
about the flexible body. Superimposing the linear 
combination of assumed modes onto the body reference frame, 
the total element displacement field can be defined. Using 
variational methods, the equations of motion of a flexible 
component in a system were written in terms of an uncoupled 
set of gross rigid and finite modal degrees of freedom. It 
has been shown (McGowan and Housner, 1985) that the 
superimpos it ion of linear elastic model onto the non-linearly 
behaving reference frame can lead to erroneous predictions, 
when either the body is highly flexible or the rotational 
speed is high. Another limitation lies in the selection of 
appropriate mode shapes which have to be consistent with the 
system configuration. This selection is further complicated 
by the configuration changes of the dynamic systems. 

In this paper, a new formulation for the transient 
analysis of multibody structural systems is proposed. This 
formulation, which overcomes the limitations stated above, is 
well suited to large space structures which consist of highly 
flexible, light weight beams. In many space structural 
applications, component beams may be connected together using 
one joint body with several hinge connections. In Fig. 1, a 


2 



typical joint body connection is shown, wherein the joint 
body JB comprises several hinges; through each of them, an 
elastic beam is connected. By adopting the finite element 
method, the total displacement field of each element can be 
defined by the three translational and three rotational 
displacements of the joint bodies to which the elements are 
hinged and the associated hinge degrees of freedom. Using 
the principle of virtual work (PVW) , the element equations of 
motion can be derived in terms of joint body displacements 
and hinge degrees of freedom. The offset of the hinge 
connection to the origin of the joint body reference is taken 
into consideration. Instead of using the constraint 
equations to govern the joint behavior, the hinge degrees of 
freedom are explicitly embodied in the formulation. 

Therefore, the kinematic inconsistency (Baumgarte, 1972), 
encountered in numerical integration of a constrained system 
of differential equations can be avoided. 

In section two, all the reference frames used to define 
the element and the joint body kinematics are introduced. 
Using the kinematic relations among different reference 
frames, the element nodal deformation with respect to the 
element corotational reference axes and the displacements of 
the element corotational axes can be written in terms of the 
displacements of the hinge connection points as well as hinge 
degrees of freedom. This is shown in section three. Using 
the virtual work, in section four, the generalized element 
internal forces which are associated with the grid point 


3 



displacements and relative hinge degrees of freedom at two 
ends are derived. In section five, the inertia matrix and 
the gyroscopic terms which include the centrifugal and 
Coriolis forces are then derived. Since there may be offsets 
from the hinge connections to the joint body reference frame, 
compensation procedures are taken in section six to rewrite 
the element equations in terms of joint body displacements. 
Accordingly, the equations of motion of different elements 
connected to the same joint body share the same generalized 
coordinates, namely the translational and rotational 
displacements of the body. 
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2. Reference Frames 


To sufficiently depict the kinematic relations existing 
among joint bodies, beam elements, and hinge connections 
between them, in this paper, six reference coordinate systems 
are introduced. These six reference frames are defined as 
follows . 



Sj, which is a Cartesian coordinate system fixed to the joint 

body, is defined such that its axes are parallel to the 

global axes in the initial configuration. As it is shown in 
Fig. 2, dj and T_. represent, respectively, the displacement 


vector of the origin of the system and the 3X3 transformation 
matrix defined in the inertial frame. Therefore, the global 
displacement of an arbitrary point p- 1 in the joint body can 
be written as 


J , , m j j 
d p = dj + Tj r p - r F 


where r^ is the position of p^ in S, 
P ^ D 


( 1 ) 


G rid point reference system (S^) 

As it is shown in Fig. 2, is fixed to a grid point p^ on 

the joint body j at which an elastic member is connected 
through a hinge. The displacements of this system represent 
the motion of the grid point. This reference is defined such 
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that its axes are parallel to which of the associated S.. 
system. Therefore, one may write 


dg = dj + Tj r«j - rg 

and 

Tg = T j 


( 2 ) 

(3) 


where d^ and represent, respectively, 
vector and the orientation matrix of S^, 
vector which defines the position of the 


the displacement 
and r- 5 is the offset 

g 

grid point in the Sj 


system . 



S -L , which is also fixed to the grid point p^ (Fig. 2), 
gh g 

defines the orientation of the hinge by using hinge axis to 

be its z-axis. Therefore this reference system is rigidly 

attached to the grid point reference system S^, and its 

y 

global orientation can be represented by 


T 


j 

qli 



(4) 


where r2 is the constant 3X3 matrix which describes the 
g h 

orientation of the grid hinge system S^ h with respect to the 

grid point reference system . 

y 

0 

F.lement corotational sys tem ( S ) 
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S^, whose motion represents the gross rigid motion of the 

element e, is defined such that its x-axis passes through two 

end nodes of the element, and y- and z-axes coincide with the 

principal axes of the cross-section of the first end node 

when the element is undeformed and twist with this cross- 

section during the dynamic motion. Standing on this 

coordinate system, the element's elastic behavior is 

0 0 

observed. As it is shown in Fig. 3, d and T are the 

c c 

translational displacement vector of the origin and the 
orientation matrix of this system, respectively. 



As it is shown in Fig. 3, is fixed to the ith (i = 1 or 2) 

node of the element e. The translational and rotational 

displacements of this system represents the motion of the ith 

end node. This system is defined such that the x-axis is 

along the tangential direction of the beam axis while the 

other two axes coincide with the principal axes of the cross- 

section of the element's end. Due to the definitions of the 
0 i 

two systems S c and S^, one may relate these two systems by 



, and due to the small nodal angular displacements resulting 
from the elastic deformation in the element, yields 
i e 1 e i 

T e = T c C T e = T c R ! i = 1, 2 ( g ) 
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where d* and are the displacement vector and orientation 

matrix of system S 1 , T 1 is the transformation matrix 

e c e 

representing its orientation in the element corotational 
0 

system S^, and 


1 i 

i -'M, 


i 

R = 


where [ 
rotations 


i i 

♦, 1 

i i 

1 


1 , <)> , < j) 1 ] 

x Y y T z 

measured in 


is 

the 


(7) 

the vector of the small nodal 

e 

corotational system S . 

c 


Elem ent Hinge system ( S A ) 

S \ . which is fixed at the ith end of the element e at which 
eh' 

the element is connected to a joint body by a hinge, is 
defined such that the hinge axis is its z-axis . Since is 

also fixed at the node, there exists a constant relation 
between this system and its associated nodal fixed system 

S . As it is shown in Fig. 3, this relation can be written 
eh 

as 


T eh ~ T e el" h (8) 

where T ^ represents the orientation of the element hinge 

system at the ith end node with respect to the inertial 

frame, and r* is the constant matrix describing the 
’ e h 
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kinematic relation between these two systems. The axes of 
the system coincide with which of the grid hinge system 

when the hinge displacement is zero. As it is shown in 

Fig. 4, resulting from the hinge displacement during the 
dynamic process, the kinematic relation between these two 
system can be written as 


i 

f eh 



v 

- 1 h 


(9) 


where represents the orientation of the element hinge 

system at ith end node with respect to the adjacent grid 
hinge system S . Since rotates about its z-axis, one may 

write 


j 


T 


i 

h 


i j i j 

cos9 h -sin9 h 0 

ij ij 

sin9 h cos9 h 0 


0 « 0 1 


( 10 ) 


where 9^ is the counter-clockwise hinge displacement from 
the initial configuration. 


9 



3. Element Kinematics 


As it is shown in Fig. 5, il and i2 are the two ends of 

element e which are connected to the joint body jl and j2 by 

hinges at grid points p-’ 1 and p-^ 2 , respectively. Standing on 

(J y 

0 

the element corotational reference frame S c , one may observe 

the small displacements at two ends due to elastic 

deformation by measuring the motion of the element nodal 
fixed systems and S^ 2 . Let [OT^ , 1 ] T and 

and [a^ 2 ,aj 2 ,a^ 2 ] T and be the 

vectors of translational and rotational displacements of the 
two element nodal fixed reference systems with respect to the 
element corotational reference axes. Because of the 
definitions in the preceding section f one may express the 
constant kinematic relations between the element corotational 
axes and the element nodal nodal fixed system as 


il 


il 


i 1 ' i? 


i 2 


il 


a v = <*„ = ct = a = oc =<k =0 


(ID 


Since the length of the element changes with time, one 
may write 

12 . 12 j2 11 11 h 

a x =1 (r 0 + d q - r D - d g ) 0 - L G (12 ) 

i 2 

where ||.|| denotes the length of a vector, r^ and r Q are the 
initial positions of the origins of the two grid point 
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where superscripts in () refer to the number of the column 
matrix in the associated transformation matrix; and, 
according to Eqs.4, 8 and 9, one may obtain 
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i 

T e 


i(l) i<2> i ( 3) 
t e t e t e 


— t r 

1 q <r h 


\ i 13 1 

t h ( 0 h ) r 


; i = il, i2; j = jl, j2 

(17) 


Therefore , in Egs. 13-16, the small nodal displacements of the 

element with respect to the element corotational reference 

axes are written in terms of the grid point translational 
displacement vectors d^ and d^ , grid point transformation 

matrices T ^ 1 and T-’ 2 , and hinge displacement B?;- 5 . 

g g n 

Since the origin of the corotational axes are fixed to 
the grid point p-* 1 , one may write the displacement vector of 

the corotational axes as 


d 


G 

C 



(18) 


And, because of Eq.6, one may write the transformation matrix 
of the corotational axes T® = [ t^ (1) , t^ (2> , t^ <3> ] in terms 

C C t— ■ 


of the transformation matrix of the first element nodal fixed 
pil 

system S as 


e Cl) _ -hU) 
t c t e 


11 —11(2) 11 — 1U3) 

<j> z t e + <!> y t e 


(19) 


e<2) 

t 


— 11 ( 2 ) 11 — 11 ( 1 ) 
t e 0 z t e 


( 20 ) 


e ( 3) _ -1K3) 

t c “ t e 


il —11 (1) 

t e 


( 21 ) 


12 


_ . . _ ^11(1) .11(2) , .11(3) 

Again, as shown in Eq.17, t , t , and t e are 

"i 1 i 2 

functions of grid point transformation matrices T J and T J , 

g g 

and hinge displacement 0^ . By Eqs . 18-21, the translational 

displacements and the orientation of the element corotational 
axes can be written in terms the motion of the grid point 
reference system as well as the hinge displacement. 
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4 . Internal Forces 


As it was defined in the preceding section, for each 

beam element, there is an element corotational reference 
frame, namely S^, with respect to which a unique elastic 

displacement field of the element can be defined. Therefore, 
in this section, the element internal forces are derived 
using the increment of the strain energy of a deformed 
element . 

Strain energy due to bendin g 

In finite element formulation, the shape functions which 
provide the displacements along the length of the beam 
element in terms of displacements at the element's end nodes 
are established. As it was used in (Housner and ef al, 1988), 
the displacement shapes for flexural motion of element e, 
whose two ends are il and 12, can be written with respect to 
the element corotational axes as 



» 

il 

. 12 


il T 12 

u c = -y c L <- 

3 il 


+ a i? <J) Z 

+ z c L r 

an <|> y + a 12 

r 

il 

12~ 



V c = L e j an § 7 

+ 

a 12 <t> z 


(23) 



il 

i2 



w c = - L e 

an 


+ a 12 <J> y 


(24) 
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where u, v, w represent the displacements along x e ~, Y e ~r and 

Q 

z -directions, respectively; x , y , and z are dimensionless 
e > r j c f j c > c 

coordinates with respect to the element length L e along the 

axes parallel to the element corotational axes, () ' denotes 
first derivative with respect to x , and 


an = (*c “ 2xg + x 3) Le (25) 

and 

a i 2 = <-x§ + xg) L e (26) 

Using the Euler-Bernoulli beam, the strains due to flexural 
motion are 


E yy = 


f f f f 

e„, = e„„ = e..„ = e. =0 


xy 


yz 


(27) 


and 


E = 


L e L 


il „ 12 

an <J> Z + a 12 0 


] + Z / L f [ 


il 


12 


an 4> y + a i2 <|) y 


(28) 


By using the constitutive equation, i.e. a = E E, one may 
obtain the increment of the element strain energy due to 
flexural motion as 



il i2 il 

5<f> z ,& |) z ,&|) y 




12 


Vr 1 [c3 


(29) 
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where 


M ’"y. * 





— 


il 

ft 

E i zA 1 1 

If 

E I 7 A 1 2 

0 

0 


4>z 

EI z a" 2 

II 

El 7 A ?? 

0 

0 


12 

<1> Z 

0 

0 

It 

E I yA 1 1 

ElyA," 


il 

*y 

0 

0 

II 

El y Al2 

II 

Ely Ap2 


i2 

<5 

L Yy J 


(30) 


in which [c*] is a time invariant matrix and 


- 


A 12 _ 


&22 ~ 


f 

1' 

1' 


1 il UA c 


a n a 1 2 dx c ~ 2 L e 


a 1? a ^2 dx c 4 Lg 


(31) 


Using Eqs.14 and 15, one may obtain the generalized internal 
forces due to flexural motion of the element as 



8u c = 8p 


(32) 


where [B^j is defined such that 


5<E> yz = L B eJ SP 6 


( 33 ) 


and 


8p = 8d 


[ 


, ji 


jl il.jl 12 j2 12, j2 

8n g 89 h 8d g 8jc g 80 h 


j2 


16 


in 

the 


which and 5 k^ are the vectors of virtual rotations of 

g g 


two grid point reference systems 


Strain energy due to stretching 

Using Eq . 13 , the mean value of the axial strain in the 
beam element can be obtained, which may be written as 


m 


e 


XX 



(34) 


As it was used in (Housner and etol, 1984), the local strain 
at the neutral axis of the element may be obtained by adding 
the first-order neutral axis strain due to flexure, which can 
be written as 


S 


E 


XX 




f 


(v c ) 


2 


+ 



(35) 


in which () denotes the first derivative with respect to x . 

c 

Therefore, the virtual work done by stretching of the element 
can be written as 


5u c - 


/. 


s s 

8e e e 

XX XX 


dv = E A L 


■/ 


S S 

8e e 

XX XX 


dx. 


(36) 


By using 


xx 




— 1 IT 


Eq.36 can be rewritten as 
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( 37 ) 


5u; 


= 8p e E ae" x l g { [b g ] + Y 2 [cf] <t> yz ] 


2 

where [c ] is a time invariant matrix 
e 


[<=:]- 


E I Z A •] i EX 7 A I? 0 0 

i » 

EI^!? EI Z A ? 2 0 0 

f i 

0 0 El yA ll El yA 1 2 

I T 

0 0 El yA 12 El y A 2 2 


in which 

An - 

Al2 = 

a 22 = 


(' ■ ■ 

I a il a il 

J 0 

I a n a i 

/ 0 

J a i? a i 
0 


dx r 


12 dx c 


15 

= - J- 

30 

= -2- L e 2 
15 


( 38 ) 


and [B ] is defined such that 
e 


r 2i i 

5a x = [b p J 8p 


( 39 ) 


Strain energy due to twisting 

Because of the assumption of the linear twisting along 
the element axis, one may write the increment of the strain 
energy due to element twisting as 
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( 40 ) 


5Uc 


GJ 

L e 


12 


^x 


i2 

<t>x 


Using Eq.16, Eq.40 can be rewritten as 




in which [B^] is defined such that 


6$, 


i? 


[b3 8p f 


(41) 


(42) 


According to Eqs.32, 37 and 41, the virtual work done by 
the total generalized internal forces which include the 
effects of element bending, stretching, and twisting can be 
written as 


5w l = 8ut 


8Uc + 8Uc 


+ 



Fi 


(43) 


The total element internal forces F? are thus obtained. 
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5 . 


Inertia 


In the preceding section the generalized internal forces 
due to the deformation of an element e are derived. In this 
section the derivation of the inertia of the deformable 
element resulting from the element large motion as well as 
elastic displacements is presented. Combining the flexural, 
stretching, and twisting effects, the total elastic 
displacement field of an element with respect to the element 
corotational axes can be written as 

ug = Nq \\f% (44) 


where 


N c 


z C L e a 1 1 - y C L G a i i z C L e a x 2 y cL e a i 1 0 

0 L e a n 0 i 2 ” Z C L e 

-L e a ix 0 ~L e a x 2 0 YcLg 


Xc 

0 

0 


(45) 


and 


T 

G 


V c - 



il 12 12 12 12 

<t> z , 4> y , <1> Z » ’ «x 1 


(46) 


Therefore, the total displacement of an arbitrary point in 
element e can be written as 


e 



d c + L g 


G 0 

(T c - T c ) 


r c + T c N c W 


(47) 
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where ()° denotes the initial conf igurat ion, and 

T 

r c =[x c ,y c , z c ). using Eq.18, and Eqs . 4 , 6, 8, and 9; 
respectively, d^ and can be written as 


d r = d 


ji 


and 


® 11 „ jl jl il _ il il 

Tc — Tg gf^ Th gT ^ R 


Let 


L e r c + N c \)/ c 


11 1* , il j- 1 

Sc = “ Tg g r h ^T h gT h 


Sc = 


ji n a jit ^ 1 ii? 

t r h r 

9 9 h il, jl 

ae u 


Sc = 


. • 1 -X 2 11 

rr il 5 jlT h ^ 

i n n* v ^ 


•■3 h 


il 

h 


ae. 


il, jl 


and 


Tc = 


S c R 


ii 


fT II 

Tc = S r R 


il 


, the virtual displacements of any arbitrary point 
e can be written as 


( 48 ) 


(49) 


on element 
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in which 


i 1 11 i l 1 

0 C = f 0 - 0 y - 0 Z 1 


Similarly, 


u g “ [ B c] P ( 51 ) 


Taking time derivative again, yields 


u g = \bI\ p + b® ( 52 ) 


where 


~ jl 

~ ji 

e 

* il 

, jl 

» * 

• 

il 

w g 

{ 

T c r 

C + 9 h 


c r c + 

S c 

R r 

• 

11, ji 

~ jl 

f * 

il 

, j 1 « 

* 

~il 

i • 

+ 0 

h 

< “g 

T c r c 

,c 

CD 

+ 

Tc 

r c + S 

C 0C 


_il 

il.il * 

~ jl 


. il, jl 

~ii 

• 

il * 

+ { 

Sc ' 

0c 

+ CO g 

s c + 

0 h 

} 0c R 

r c 
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- jl 

+ 2 { CO g Tc N 


.e .11, jl 

Vc + 0h T < 


• e 


.11 


11 


. e 


N c \|/ c + S c 4> c R N c V|l c } 
(53) 


Since the virtual work done by the inertial forces of 
element e may be expressed as 


5w^ = p 



(54) 


0 1 

where p is the mass density; using Eqs.50 and 52 r one may 
write 


5w l = 5 P e [M e ] 


. . e 

p + 



(55) 


where [M e ] is the element inertia matrix and g e is a vector 
of quadratic velocity terms which include the centrifugal and 
Coriolis forces of' the element, which can be written, 
respectively, as 


and 


e e 

[M ] = p 


e ° 

g = p L e 


L e 




G G 

[B c ] b c 


[Be] 


dx c 


dx c 


(56) 


(57) 
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Using Eq.43 and 55, the total virtual work of the 


element can be written in matrix form as 

X e 

5w t = 5p {[M]p + g + F i - Q } = 0 


(58) 


where Q e is the vector of generalized external forces which 
can be obtained by using the virtual work. 
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6. Offset Compensation and Connectivity 


In the previous sections, the element equations of 


motion are derived in the generalized space spanned by the 


displacements (translational 
point reference systems 


and rotational) of the grid 

and S-^ on the joint bodies where 
9 


the element is connected, as well as the hinge degrees of 
freedom existing between joint bodies and the element. In 
many cases, a joint body may be used to connect several 
elements. Using the displacements of each connecting grids 
on the same joint body to be the system's generalized 
coordinates, the efficiency in solving the dynamic equations 
must be decreased due to the increased number of dependent 
generalized coordinates. Therefore, in this section, a so- 
called offset compensation method is introduced. On each 
joint body, only one grid point, which will be referred to as 
the primary grid point, is used to represent the motion of 
the joint body. All the other grid points on the joint body 
where elements are connected are referred to as the offset 
points to the primary grid point . Using the PVW, the 
equations of motion of each element connected to the offset 
points can be rewritten in terms of the motion of the primary 
grid point as well as the hinge displacements. The system's 
equations of motion can be thus assembled by using the 
minimum number of generalized coordinates. 
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Using the kinematic relations between an offset point 
and its associated primary grid, which is as written in Eqs . 2 


and 3, it yields 

~ e e _ e 

5p = [L ] 8q 


(59) 


where 

T 

8q S 

and 


e 

[L ] 


8dji 8 ji . 

11. jl 

Sts* 

8d j2 

12,j2 1 

Sitjj 89 h J 

I 3 

- (Tjir’ 1 ) 

0 


0 

I 3 

0 

0 6 

0 

0 

1 

I 3 - (T j2 r g 2 ) 


0 6 


o I 3 


0 0 


0 

0 

1 


in which O^'s are 6X6 null matrices, and are 3X3 identity 
matrices. Similarly, 


. e 

p 


[L ] q 


(60) 


Taking derivative with respect to time, one may get 


..e e ,.e e 

p = [l ] q + x (q,q) (6i) 
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where 


ji 

- 0) CO (T jir g ) 

0 

0 

X - ~3* _ i2 - 

- CO CO (T j2 r g ) 

0 

0 

, in which CO's are vectors of angular velocities. 

Substituting Eqs.59 and 61 into Eq.58 f the element 
equations can be written as 


[M ] q = F (62) 


in which 


— — s 

[M ] = [L S ] [M 6 ] [L e ] 

and 

[F 6 ] = [L S ] { Q® - Ft - g 6 - [M 6 ] %} 
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Since all the elements connected to a joint body share the 
same generalized coordinates, which are the translational and 
rotational displacements of the joint body, the system 
equations of motion can be assembled by the element equations 
of Eq.62 using the minimum number of generalized coordinates 
with minimum number of constraint equations. 
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7 . Summary 


In this paper, a finite element approach is used to 
formulate the dynamic equations of joint-dominated flexible 
multibody systems undergoing large rotations. On each 
element, a set of element corotational axes are selected to 
represent the rigid modes of the element motion. Using the 
element kinematics defined in the paper, the element rigid 
modes and the small elastic deformation with respect to the 
element corotational reference can be written in terms of the 
displacements of the grid points on the joint bodies, at 
which the two ends of the element are connected through 
hinges, and the associated hinge rotations. Using the PVW 
(principle of virtual work) the element equations of motion 
can be derived. In this method, the hinge connection is 
defined by introducing a relative degree of freedom instead 
of using the constraint equations which may cause the system 
configuration inconsistency due to the truncation errors 
accompanied with the numerical integrations. To further 
reduce the number of generalized coordinates in the system, 
the offset compensation method is used for the element 
connectivity at the joint body. Based on the formulations 
presented in the paper, a general purpose transient dynamics 
code LATDYN (Large Angle Transient Dynamics) has been 
developed in NASA Langley and used for the design and 
analysis of space crafts and structures. Even though only 
beam elements and hinge joints are used presently, the method 
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presented in this paper may also be applied to more general 
systems consisting of different elements as well as various 
joint connections. This is worth further studies. 
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Figure 2 : Reference Frames on Joint Body 
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Figure 3 


Reference frames on a Beam Element 


34 



i 

z 



gh 


Figure 4 : Hinge Degree-of-Freedom 
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